Unrestricted renormalized mean field theory of strongly correlated electron systems 
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We generalized systematically the renormalized mean field theory in the case of uniform states to 
the unrestricted case of general inhomogeneous states with competing spin-, charge- and supercon- 
ducting orders. Applying the theory to high-T c superconductors, we discuss the issues of electronic 
inhomogeneity, the superfluid density, and in particular the local electron density of states. The 
results account for many intriguing aspects of the phenomenology. 



in 



a ■ 
o ■ 

■ 

1 1 

5h ■ 

Or 

s : 

. 



i 

-a 
a 
o 

o 



> 
(N 

o 

O 

in 
o 

1 



i 

-a 
c 

o 
o 



As a prototype of strongly correlated electron sys- 
tems, high-T c superconductors reveal rich phenomenolo- 
gies. Apart from the apparent competing anti- 
ferromagnetic and superconducting orders, recent scan- 
ning tunnelling microscopy (STM) measurements have 
suggested checkerboard charge-density-wave order in 
under-doped NaCaCuOCl and BiSrCaCuO superconduc- 
tors 0]. On the other hand, earlier STM measurements 
already revealed nano-scale spatial inhomogeneity of the 
quasi-particle energy gap Q, but yet the local density of 
states (LDOS) at very low energy are remarkably homo- 
geneous, seemingly ignorant of the gap inhomogeneity Q. 
The physics behind the robust low energy quasi-particle 
states are intriguing. In short, strongly correlated elec- 
trons may support complex phases and their response to 
extrinsic impurities may be highly nontrivial. This moti- 
vates us to develop a reliable effective theory that takes 
proper care of the strong correlations and allows compet- 
ing orders and inhomogeneous electronic states from the 
starting point. 

A widely adopted model for high-T c superconductors 
is the t-J model, 



t-j 



h.c. 



J5> 



Si 



(1) 



where the constraint of no double occupation at any site 
is implied. In order to fit the band structure, hopping 
terms connecting next-nearest-neighbor sites and further 
neighboring sites may be necessary but they do not af- 
fect the following development of the theory. Recently 
the t — U — J model has also been used to study the Gos- 
samer superconductivity 0, in which the constraint of 
no-double-occupation is relaxed but a Hubbard [/-term 
is added, 



H = H t -j + [/> y n^n^;. 



(2) 



The following development is formulated within this gen- 
eral model, Eq.©, with obvious specific limits to the t— J 
and Hubbard models at our disposal. 

A very powerful way to take proper care of the strong 
correlation effects is to Gutzwiller project a trial wave- 
function, = P\ip ), where \ifto) denotes a free-particle 



many-body wave-function, with possible order parame- 
ters, and P = Hi(l — aDi) is a projection operator. 
Here Di = n^riii, and a is a variational parameter be- 
tween and 1. The calculation can be done by Monte 
Carlo methods |5(. However, it is limited by the lattice 
size and the parameter space. Alternatively the effect 
of Gutzwiller projection is taken into account approxi- 
mately in a re-normalized mean field theory (RMFT) |(| , 
which agrees qualitatively with the Monte Carlo re- 
sult. The Gutzwiller projected d-wave BCS states ap- 
pear to capture the basic features of high-T c phase di- 
agram 0, 0, 0- I n this paper, we develop a system- 
atic extension of the RMFT to the case of unrestricted 
spin/charge densities. We then discuss the physical quan- 
tities based on the unrestricted RMFT and explain the 
robust low energy states observed in the STM. 

Let |^o) be a Hatree-Fock state which carries the bare 
variational unrestricted order parameters. In the usual 
mean field theory this wave function is used to calcu- 
late the expectation value of the Hamiltonian, and the 
optimization with respect to the order parameters leads 
to the mean field self-consistent equations. For strongly- 
correlated electrons, however, this usual type of approach 
misses the most important physics, namely, the possi- 
ble emergence of Mott insulating phases. This physics 
is conveniently picked up by Gutzwiller projecting the 
Hatree-Fock state, namely, the trial state is given by 
|V-') = P |V^o)- However, we are no longer justified to use 
a projector P with identical projection parameter a at 
each site, since the effect of strong correlations depends 
on the electron density. This motivates us to use a gener- 
alized projection operator V — HiVi so that = V\ipo), 
with 

Vi = y?(l - on Di) = Et + yiQi + rm%D u (3) 

where E l = (\ ~ n iT )(l - rc.jjj and Q, = J^a- n w{^ - n ia) 
are the empty and single occupation operators, respec- 
tively. These are standard projection operators, satisfy- 
ing E? = Ei, Qf = Qi and D\ = D L . Therefore 



(4) 



Here rji = 1 — a,, and the fugacity fji is introduced so 
that the local charge density is not changed by projec- 
tion iHEii or fi = (V>ok|^o> = (^klVVWW, the 
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advantage of which is that the variation of rji does not 
influence the local charge density. (In the t — J model 
rji = 0.) We note that we are working in the grand- 
canonical ensemble, which proves more convenient than 
the usual canonical ensemble approach usually used in 
the case of uniform state RMFT. 

We shall use the shorthand notations (-) = (i/'ol ■ \i>o) 
and (•) = Z (tpo\V ■ V\ipo) with the normalization fac- 
tor Z = {iJjo\V 2 \4>o} = (P 2 )o- Since \ip ) describes free 
particles, one can use Wick's theorem to obtain all the 
contractions in Z . Under the Gutzwiller approximation, 
inter-site correlations are ignored when expectation val- 
ues of projection operators are evaluated. One thus ob- 
tains Z = TliZi, with 

z i = CPi)o = eoi + yfqoi + ViVidoi, (5) 

where e 0i = (Ei) Q = 1 - /, + d Ql , q oi = (Qi) = f. L - 2d oi 
and doi — (Di)q = r^r^. In the above derivation, we 
have assumed the absence of on-site pairing and spin flip- 
ping order parameters in the wave- function IV'o)- Simi- 
larly, we have a = = e i/zu q t = (Qi) = qoiyf/zi, 
and di = (A) = dmrjfyf/zi. The fugacity y l is de- 
termined by enforcing /, = qi + 2dj, and is a function 
of rji, fi and d^ (or equivalently tv). The calcula- 
tion of yt is even unnecessary. Eliminating we find 
diei/qf = rfdoieoi/qQi, which determines di uniquely. 

We now discuss the effect of projection on the local spin 
moment directed in the z-direction. Since Si is already 
projective, we have ViSiPi = yfSi, and therefore 

m i = (Sf) = (S^oyf/zi = g 8 (i)m 0i , (6) 

where g s (i) = qt/qoi is the re- normalization factor for 
the variational spin moment moi = (S£)q = (r^ — t\l)/2. 
Similar consideration applies to the spin-spin exchange, 

(Si ■ Sj) = g s (i)g s (j)(S l ■ Sj) . (7) 

Using the identities 

ViCujVi = [yi{l - n lS ) + r)iVinicr]Cia, (8) 
ViClVi = Ml - n lS ) + myfn m ]Cl, (9) 

we obtain 

(ClC ja ) = g t a(i)gta(j)(Clc jtT ) , (10) 
with the re-normalization factors 

9t<r(i) = (1 + n*\ -^j-- (ii) 

V eoiQoi y qoidoi 

In Eg. 1)11(1 the first term arises from the hopping process 
that does not encounter a-electrons at site i, and the sec- 
ond is from that involving a double occupation. We note 
that the re-normalization factors for spin-spin exchange 
and the hopping reduce in disguise to the known results 



in uniform charge density and uniform magnetically or- 
dered states Q . The result in the non-magnetic case for 
the t — J model was used for granted in Ref. 0] . 

The total internal energy E — (H) can then be written 
as, 

E = -t 9ta{i)gta{j){x* ll( j + Xijv) 

3 J 

(v) « 

where we have defined Xija — (C\ a Cj a )o and Xij — 
T, a Xija, Ay = (Bij)o with Bij = C^C^ - C^C n , and 
used the standard decomposition 

3 

(Si ■ Sj) = m 0l m 0j - -(xtjXij + A* 3 Ay). (12) 

It is understood that r ia , Xija an d Ay are determined 
by tpo, di is determined by and r]i, and therefore E is 
eventually a functional of ipo and {rji}. 

The optimization of E requires SE/S(tf>o\ = and 
dE/drji = 0. The first condition leads to a Schordinger 
equation Hmf\^o) — A|^o)) with the mean field Hamil- 
tonian 

+ E(^+^)+E(£~m)^, (13) 

(y) «<T 

where \i is the chemical potential. This is a self-consistent 
unrestricted RMFT in the sense that 1) ^o) is deter- 
mined by Hmf, and so are the variational order param- 
eters; 2) There is no restriction on the spatial variation 
of the order parameters; 3) The effect of strong corre- 
lation is reflected in the re-normalization coefficients gt a 
and g s ; 4) The re-normalization coefficients depend on rji, 
and the global minimum of the energy E is reached af- 
ter the self-consistency of the mean field theory and the 
optimization of rji are simultaneously achieved. In the 
case of t — J model one simply sets rji — ,12|. The the- 
ory has the advantage that the Mott physics is built in. 
In the literature ^j^i the unrestricted Gutzwiller projec- 
tion has been applied to the multi-band Hubbard models 
in dealing with para- and ferromagnetic states in transi- 
tion metals. Further encoded in the present theory are 
the intriguing unrestricted competing anti-ferromagnetic 
and superconducting orders, which are most interesting 
in layered cuprate superconductors. 

The mean field internal energy (Hmf + M-^V)o is dif- 
ferent from E, but a self-consistent minimizes E so 
that SE = J2n £ nSf n + J2 n .n> K,n'^/n<^/n' H where 
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{e„} is the single-particle spectrum of Hmf, fn is the 
Fermi-Dirac occupancy of the n-th single-particle orbital, 
and V n>n i is the residual two-body interaction kernel. 
The mean-field single-particles are not necessarily the 
electrons themselves, but are the bare Landau quasi- 
particles. The interpretation of these states, and in par- 
ticular, the excitation energy spectrum, would be ob- 
scured if the charge density in \iJjq) were different from 
that in \ip). 

In what follows we shall apply the theory to the t — J 
model, relevant to the cuprates. First, we discuss the ef- 
fect of the projection on the observables. Since the the- 
ory is projective, the order parameters in the RMFT are 
different from the measured ones. The re-normalization 
in the spin moment nij = g s (i)moi has been derived 
previously. For the superconducting pairing order pa- 
rameters, it can be shown easily using the method de- 
scribed above that A» 3 - = (-By) = g&(i, with 
9A(i,j) = \ Y, a 9ta{i)gta{i)- In the charge-uniform non- 
magnetic states, 5a = 2x/(l + x) < 1 where x (with — 
x) is the hole-doping level away from half-filling. This 
result is well-known in the literature, and is argued to be 
related to the existence of pseudo-gap 0, Q • The pair- 
pair correlation function on disconnected bonds is re- 
normalized in a similar fashion, (B^Bki) = g\{B\^Bki)o- 
This justifies the assignment of A as representing the 
long-range off-diagonal order. However, the pair-pair 
correlation on the same bond is re-normalized quite dif- 
ferently, since B\^Bij projects out a singlet spin pair 

on the bond (ij), (B^Bij) = —2(Si ■ Sj — jriiTij) = 

-2g s (i)g s {j){S i -S :j )o + \f l f r (In the t—J model = n t 
so that qi = fi.) Therefore we expect significant incoher- 
ent part in the pair-pair susceptibility due to short-range 
spin correlations. 

The superfluid density p s is related to the second or- 
der response of free energy to a vector potential. At zero 
temperature, it is related to the kinetic energy and is 
therefore re-normalized in the same way as the hopping 
parameters are, so that p s — gl a tx — 2xt\/ '(1 + x) in 
a uniform state, where \ is the variational hopping or- 
der parameter in the RMFT. (The expression is slightly 
changed in the presence of longer-range hopping, but the 
re- normalization factor is the same.) One realizes that 
this relation is consistent with the result of elaborate 
gauge theory 0|. Indeed integrating over gauge fluc- 
tuations in the gauge theory effectively restores, to some 
extent, the strong correlation effects. The appealing fea- 
ture of the projected density functional theory is that no 
slave degrees of freedom, and therefore no gauge fields 
are involved. 

The factor gta{i) measures the overlap between the 
bare quasi-particle state in the RMFT and a corre- 
sponding real electron wave function at the given site. 
Therefore, the electron Green's function may be written 
as G(ia,ja) = g t <j(i)gtati)Go( i<7 , 3 a ) + G inc {ia,ja) 0, 



where the first contribution is the coherent part from 
RMFT. We suppressed the energy arguments for brevity. 
Consider the on-site Green's function. A corresponding 
decomposition applies for the spectral function, A a (i) = 
g^(i)A 0a (i) + A inc>a (i). Even though RMFT does not 
provide Ai nCtCr directly, some properties of it can be ar- 
gued, as discussed previously in the literature 0|. We 
extend the discussion here for generally inhomogeneous 

cases. Since /^Ea^W = f-oc £<r A o*(i) = fi 
but gta{i) < 1, one concludes that there is incoherent 
spectral weight in the occupied side. While the gen- 
eral expression is quite complicated in general, it has a 
simple form for a nonmagnetic site, which is given by 

/"ooEo-^cW = (1 - e 07(l + ei). On the other 
hand, the number of unoccupied states at site i is given 
by Eo- J duj{Ci a 5{ijj — H)C\ a ) — 2e i7 where we dropped 
a contribution qi from the upper Hubbard band , which 
is out of the Hilbert space of the t — J model. To this 
total unoccupied spectral weight, the contribution from 
the coherent part of G, which does describe excitations 
below the Mott gap, is given by 

?^ )2(1 - riCTH l + 4^1-e?)- 

Of particular interest is the nonmagnetic case (moi = 0) , 
in which the above discussion predicts that the unoccu- 
pied spectral weight is exhausted by the coherent part, 
in agreement with a recent argument 0, . By conti- 
nuity, one would expect that the incoherent part in the 
occupied side vanishes while approaching the Fermi level. 
The fact that incoherent excitations are abundant in the 
occupied side preferably at higher energies implies that 
1) the nodal excitations at the Fermi energy are quite 
robust against strong correlations and 2) the antinodal 
excitations are largely incoherent due to strong correla- 
tions. This seems to be consistent with the the energy- 
dependence of the spectral function widely observed in 
angle-resolved photo-emission measurements |l8| . 

We now calculate the LDOS in the t—J model with 
weak random on-site potential, relevant to STM data. 
The inset of Fig. I (a) shows the scalar impurity potential 
profiles. The main panel of Fig. 1(a) shows the coherent 
part of the LDOS calculated from the RMFT at an av- 
erage doping level x — 0.1. The system carries d-wave 
pairing order without spin moments. While the DOS 
at and above the gap energy scale disperse from site to 
site, the low energy DOS seems to be independent of 
position and ignorant of the random impurities, in qual- 
itative agreement with experiments Q. In order to see 
whether this is a universal feature of doped Mott insula- 
tors, we deliberately erased the pairing order to find that 
the system evolves to a staggered-flux (or d-density-wave) 
state, with circulating currents whose chirality switches 
from plaquette to plaquette 0. The LDOS in this state 
is shown in Fig. 1(b). Now the low energy LDOS at dif- 
ferent sites no longer merge together, but those at an 
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energy above the Fermi level does. Since the DOS at 
the band-gap minimum vanishes linearly with energy in 
a charge-uniform staggered-flux state in the same way as 
at the Fermi level in a uniform c£-wave paired state, we 
conclude that the homogeneity in the low energy DOS 
in a <i-wave superconductor is merely a result of vanish- 
ing DOS in the parent uniform phase. In other words, 
the strong correlation effect stabilizes the d-wave pairing 
or density-wave states, and then some quasi-particles are 
protected by the small phase space available to elastic 
scattering. 

A few remarks are in order. First, we are exaggerating 
the density of impurities in Fig.l. The phase-space pro- 
tection of the quasi-particles is limited to weak disorder 
and low impurity concentrations. At stronger disorder 
and denser impurities, the inhomogeneities of the LDOS 
prevail at all energies. Second, in comparison with the 
U(l) slave boson mean field theory [2fJ, our results are in 
better agreement with experiments. This is because the 
U(l) mean field theory underestimates the kinetic energy 
so that the holes are more susceptible to impurity local- 
ization. In addition, as compared to calculations with- 
out the projection and the self-consistency, the LDOS in 
our case are more uniform against potential impurities 
in space. The physics is best seen in the Mott insula- 
tor limit, where the low- lying excitations are spinons, 
which are completely ignorant to potential impurities. 
This means that strong correlation effect reduces the 
charge susceptibility. Third, the peak-to-peak gapin our 
case varies less severely than in experiments |2j, |3( . We 
note that in order to mimic the experimental data, gap- 
impurities were introduced phcnomcnologically in a re- 
cent model j2jl ■ It would be highly interesting to find a 
mechanism that could lead to such gap impurities. 
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